Experimental magnetic form factors in C03V2O8: A combined study of ab initio 
calculations, magnetic Compton scattering and polarized neutron diffraction 
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We present a combination of ab initio calculations, magnetic Compton scattering and polarized 
neutron experiments, which elucidate the density distribution of unpaired electrons in the kagome 
staircase system Co3V208- Ab initio wave functions were used to calculate the spin densities in real 
and momentum space, which show good agreement with the respective experiments. It has been 
found that the spin polarized orbitals are equally distributed between the t2g and the Cg levels for 
the spine (s) Co ions, while the Cg orbitals of the cross-tie (c) Co ions only represent 30% of the 
atomic spin density. Furthermore, the results reveal that the magnetic moments of the cross-tie Co 
ions, which are significantly smaller than those of the spine Co ions in the zero-field ferromagnetic 
structure, do not saturate by applying an external magnetic field of 2 T along the easy axis a, but 
that the increasing bulk magnetization originates from induced magnetic moments on the O and 
V sites. The refined individual magnetic moments are /i(Coc) =1.54(4) /is, /i(Cos) =2.87(3) /is, 
^(V)=0.41(4) fiB, /i(Ol)=0.05(5) fiB, /i(O2)=0.35(5) /is and /i(O3)=0.36(5) /is combining to the 
same macroscopic magnetization value, which was previously only attributed to the Co ions. 

PACS numbers: 75.25.+z; 75.30.Et; 75.30.Gw 



I. INTRODUCTION 

C03V2O8 represents the 3d transition metal ortho-oxo- 
vanadates labeled as kagome staircase structures and 
crystallizes in the orthorhombic space group Cmcai*^ 
Its crystallographic structure is characterized by edge- 
sharing CoOe octahedra forming buckled layers of corner- 
sharing triangles, the kagome staircases, which are sepa- 
rated along the b axis by VO4 tetrahedra (Fig. [1]). The 
magnetic exchange is effectuated by a 90° Co-O-Co in- 
tralayer pathway. Interestingly, this system exhibits a se- 
quence of five magnetic phase transitions featuring tem- 



perature dependent incommensurate antiferromagnetic 
phases with commensurate lock-ins and a ferromagnetic 
ground state, where all magnetic structures are collinear 
along the a axis.^ The ferromagnetic structure reveals 
two strongly different magnetic moments for the cross- 
tie Co ions (4a site) and spine Co ions (8e site) of 1.54 
jj^B and 2.73 /iB,^ respectively, despite the fact that both 
Co^+ ions apparently present high-spin configurations 
as macroscopic measurements exhibit saturation of the 
cross- tie moments In the kagome staircase structure 
several nearest and next-nearest neighbor exchange in- 
teraction pathways are possible^ which is the motivation 
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FIG. 1: (Color online) Visualization of the kagome staircase 
structure of C03V2O8. (a) shows the edge sharing MOe octa- 
hedra [CocOe light blue (light gray), CosOe dark blue (dark 
gray)] isolated by non-magnetic VO4 tetrahedra (gray), (b) 
depicts a single kagome staircase viewed along the b axis with 
only the magnetic ions on both crystallographic sites [Coc 
light blue (light gray), Cos dark blue (dark gray)]. 

for this study. Investigating the magnetization density, 
which is reflected by the respective magnetic form fac- 
tors of the ions, may reveal preferred exchange pathways 
with the presence of magnetization on the involved O 
sites. Eventual induced magnetization on the empty d- 
shell of V sites would allow interlayer coupling by super- 
super exchange. Therefore, magnetic Compton scatter- 
ing and polarized neutron diffraction experiments have 
been carried out leading to the spin density in momen- 
tum space and the magnetization density in real space, 
respectively. However, detailed and precise analysis of 
these quantities are required in order to determine the 
exact contribution of the atomic species involved in the 
studied system. Quantum chemical modeling is there- 
fore needed to gain insights, at a molecular level, into 
the electronic structure of the two cobalt-oxide octahe- 
dra. In this context, ab initio cluster calculations were 
performed for the CocOe and CogOe octahedra yielding 
precise molecular orbitals (MO) and wave functions (wf). 
The latter were used to analyze the experimentally ob- 
served density distributions. The refinement of the con- 
tribution of each MO/wf, at a quantum chemical level, 
to the real and momentum space densities simultaneously 
is a powerful procedure allowing interesting and valuable 
features of the magnetic form factors to be determined. 



II. AB INITIO CALCULATIONS 

The two different clusters CocOe and CogOe were 
modelled separately. The calculations were performed 
within the framework of the Kohn-Sham formulation of 
the density functional theory using the PC GAMES S 
program.- The functional B3LYP was employed to 
approximate the exchange-correlation interaction. The 
B3LYP is a hybrid functional, well adapted to study 
transition metal compounds and magnetic interactions. 



in which a predefined amount of the exact Hartree- 
Fock exchange is added to the well known pure density 
functionals^^iii^i^ The atoms in the cluster were described 
using Ahlrich's pVDZ atomic orbital (AO) basis sel3 
Co(14s,8p,5d,lp)/[5s,2p,2d,lp],0(7s,4p,ld)/[3s,2p,ld]. 
The notations (klm) and [klm] indicate the number of 
Gaussian type orbitals and contracted Gaussian type 
orbitals, respectively. In order to mimic the Madelung 
potential, the two quantum mechanical clusters were 
surrounded by point charges (PC) according to the 
Effective Fragment Potential methodJ^ As previously 
reported for other systems ^li^ii^ii^ the choice of the 
embedding method was shown to be crucial for the 
physical meaning of the ab initio calculations. To 
avoid the electron density leaking out of the cluster, a 
boundary region has been introduced, which is formed 
by effective core potentials (ECP) placed in the nearest 
cationic positions around the cluster. Thus, the first 
coordination shell of Co^+ and V^^ ions has been 
described by ECPs according to the SBKJC ECP basis 
set.^^ As native ECPs for Co and V do not treat the 3s 
and 3p electrons as core electrons, those ECPs would 
in fact be too compact due to the higher number of 
valence electrons. In order to overcome this problem 
the ECPs for Mg and Al have been used due to the 
fact that their ionic radii are closer to those of Co and 
V, respectively. Except for the 3d shells, the remaining 
electrons are replaced by an effective potential. Thus, 
1565 PC (3x3x3 unit cells with the respective Co ion 
in the center) have been built to mimic the Madelung 
potential on the cluster. The MO coefficients relevant to 
the Co3d were extracted from the simulations and used 
to model the Magnetic Compton Profiles (MCP) and 
the magnetic form factors. 



III. EXPERIMENTAL 

A. Unpolarized neutron diffraction 

C03V2O8 single crystals have been grown from self- 
flux in a Zr02/Y crucible by the slow cooling method. 
Some of the crystals have been ground and investigated 
at the high-flux neutron powder diffractometer D20 (In- 
stitut Laue Langevin, Grenoble) confirming the absence 
of parasitic phases. The nuclear structure of a chosen 
single crystal has been studied at the four-circle diffrac- 
tometer D9 (Institut Laue Langevin, Grenoble). A set of 
more than 500 independent reflections up to sin ^/A=0.92 
has been measured using two different wavelengths in 
order to determine the extinction effects (Ai =0.835 A, 
A2=0.512 A). The data collection has been performed in 
the paramagnetic phase at T=13.5 K, which is just above 
the Neel temperature of 11.2 K. In order to reveal pos- 
sible structural changes between the paramagnetic and 
the ferromagnetic phase due to structural phase transi- 
tions or magnetostriction an additional nuclear structure 
investigation has been carried out at the single crystal 
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diffractometer D15 (Institut Laue Langevin, Grenoble) 
under the same experimental conditions as the polarized 
neutron diffraction experiment (Sec. IIII B]) . i.e. at T=3.5 
K with an applied magnetic field of H=2 T along the 
easy axis a and using a wavelength of A=0.854 A. 

B. Polarized neutron diffraction (PND) 

The real space magnetization density has been studied 
at the hot neutron spin polarized two-axes diffractometer 
5C1 (Laboratoire Leon Brillouin, Saclay). Neutrons from 
the source are monochromated and polarized by the (111) 
refiection of a magnetized Heusler crystal Cu2MnAl. The 
wavelength is 0.84 A, which corresponds to the maximum 
fiux of the hot source and is ideal for studying large do- 
mains of reciprocal space. The polarization factor of the 
beam is p = —0.88. In order to fully magnetize the sam- 
ple and avoid beam depolarization a magnetic field of 
H=2 T has been applied along the easy axis a. The 
fiipping ratios R (Eq. [1]), the ratios between the spin- up 
and spin-down intensities, of over 500 independent {hkl) 
refiections with /i=0,l,2 have been measured in the fer- 
romagnetic phase at T=3.5 K. 

^ ^ (F^ + q^Fl)p+ + 2FNq^FMpt + (1 - Q^)Q^Fl,ypm 
i^N + Q'FIi)Pp + 2Fm'FMP^ + (1 - q')q'Fi,yp^ ' 

(1) 

Fm and Fn denote the magnetic and nuclear structure 
factors. q=sma is a geometric factor with a being the 
angle between the scattering vector and the magnetiza- 
tion vector. The parameters p^^^ and i/pm are extinction 
correction factors of the respective cross-sections accord- 
ing to Ref. [13. 

C. X-ray magnetic Compton scattering (MCS) 

The investigation of the spin density in momentum 
space has been carried out at the High Energy Inelas- 
tic Scattering beamline BL08W at SPring-8 in Hyogo, 
Japan. This beamline is designed for magnetic Compton 
scattering spectroscopy as it offers high energy ellipti- 
cally polarized X-rays emitted from an Elliptic Multi- 
pole Wiggler. The incident photon beam with an energy 
range of 170-300 keV is monochromated and focused by 
an asymmetric Johann type monochromator using the Si 
(620) refiection. The sample magnetization is achieved 
with a superconducting magnet with a maximum field of 
2.5 T and a minimum polarity-switching time of 5 sec- 
onds. The backscattered photon energy is analyzed by 
a 10-segmented Ge solid state detector positioned at a 
scattering angle of 178.4°. The experiment has been car- 
ried out with an incident photon energy of 176.3 keV, 
which gives a good compromise between the beam inten- 
sity and the scattering cross section. The initial interest 
of applying this method to C03V2O8 was to map the mo- 
mentum space spin density of the ferromagnetic phase as 



a projection onto the b*-c* plane in order to gather infor- 
mation about the 3d electron spin states and to correlate 
the results with those obtained from the polarized neu- 
tron diffraction experiment, but the experimental con- 
ditions and especially the large magnetic anisotropy of 
the system did not allow that. The minimal achievable 
sample temperature is approximately 5.6 K, i.e. close 
below the magnetic transition into the antiferromagnetic 
phase. It can be seen in the magnetic phase diagrams of 
Co3V2 08^'^^ that at this temperature already weak mag- 
netic fields applied along the 6 or c axis induce a mag- 
netic phase transition into the antiferromagnetic phase, 
while H\\a stabilizes the ferromagnetic one. The neces- 
sity of applying a magnetic field of considerable strength 
and therewith magnetizing the sample along the inci- 
dent beam in order to increase the magnetic contribution 
to the scattering cross-section beside the requirement of 
turning the sample about a vertical axis to be able to 
reconstruct the two-dimensional momentum density, led 
to a change of strategy. To make sure not to induce 
magnetic phase transitions by rotating the sample with 
respect to the field direction the measurements have been 
carried out within the antiferromagnetic phase at T=7.5 
K and by applying a magnetic field of H=2 T with the 
induced ferromagnetic component lying in the b-c plane. 
In addition to the trivial directions [010] and [001] four 
further directional MCPs Jmag in the Py-Pz plane have 
been investigated. An additional profile has been mea- 
sured along the [100] direction by decreasing the applied 
magnetic field to 0.25 T. A directional MCP yields the 
projection of the spin momentum space density onto the 
scattering vector, which is by definition p^: 

JmagiPz) = J J |Xt(p)P " \Xi{p)\^ dp^dpy (2) 

— oo — oo 

with XT(i)(p) denoting the momentum wf of an occupied 
majority (minority) spin state. The Compton profiles 
of the respective sample magnetization states have been 

recorded for 60 seconds, repeating the cycle [H \ h 

H — ] multiple times. 



IV. RESULTS 

A. Nuclear structure and extinction 

The nuclear structure investigation at 13.5 K on 
the four-circle diffractometer D9 confirmed the correct 
phase formation of the orthorhombic structure (space 
group Cmca) with a=6.015(3) A, 6=11.480(5) A and 
c=8.289(4) A. The observed integrated intensities have 
been corrected for absorption by applying the transmis- 
sion factor integral exp[/i(tin H-fowt)] and analyzed by si- 
multaneously fitting a structure model to both datasets 
with Ai and A2 using FullProf.^^ The refinement process 
included the atomic positions and isotropic temperature 
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factors of Co and O plus three additional extinction pa- 
rameters according to an anisotropic Shelx-like empirical 
extinction correction*^ The atomic position and temper- 
ature factor of V have been fixed in all refinements, be- 
cause of its low coherent neutron scattering cross section. 
The refined structural parameters (i?=5.3) are listed in 
TabUl 

In order to analyze the nuclear structure under identi- 



TABLE I: Structural parameters of the investigated C03V2O8 
single crystal. 



Atom 


X 


y 


z 


^(A^) 


Col 











0.27(10) 


Co2 


0.25 


0.1328(7) 


0.25 


0.20(7) 


V 





0.3773 


0.1204 


0.30 


01 





0.2489(3) 


0.2700(4) 


0.38(4) 


02 





0.0008(4) 


0.2448(4) 


0.33(4) 


03 


0.2702(3) 


0.1185(3) 


0.9990(2) 


0.33(4) 



Extinction parameters 
xii=1.0(l) X22=0.36(5) X33=0.6(l) 



cal conditions as the fiipping ratio measurement (T=3.5 
K, H=2 T), only those refiections have been measured 
for which the ratio 7 between the magnetic and the nu- 
clear structure factor has been derived from the observed 
fiipping ratios. The ferromagnetic contribution to the in- 
tegrated intensity / has been canceled out according to 



N 



IQmI 



F? 



N 



M 



Fj,{l+qW), (3) 



where Qm is the magnetic interaction vector. The 
subsequent refinement showed that no considerable 
change in the nuclear structure has taken place, i.e. the 
derivation of the Fm from observed fiipping ratios by 
using the observed Fn at T = 13.5 K is justified. 

All low-angle nuclear refiections suffer considerably 
from extinction, therefore, special attention has been 
paid to the extinction of magnetic scattering. As the fiip- 
ping ratio treatment uses the same extinction parameters 
for both nuclear and magnetic scattering, it is important 
to verify, if the extinction effects are indeed comparable. 
Therefore, three strong magnetic refiections have been 
measured as a function of applied magnetic field after 
the sample has been cooled in zero-field to 3.5 K. Fig. [2] 
shows the integrated intensities of three refiections after 
the nuclear contribution has been subtracted. The field 
dependence of the magnetic contribution reveals a sur- 
prising and interesting tendency: Instead of increasing 
with increasing applied field, as one would expect if the 
cross-tie moments get saturated, the intensity of mag- 
netic scattering drops significantly. This observation can 
be explained with field dependent increase of primary 
extinction. At H={) T the sample exhibits a multido- 
main state with presumably negligible extinction effects. 




FIG. 2: (Color online) Intensity of three different magnetic re- 
flections in dependence of an applied magnetic field revealing 
primary extinction effects. 



By increasing the field the magnetic domains grow until 
they reach approximately the size of the structural do- 
mains. On reaching saturation at H^{).2b T the primary 
extinction effects for magnetic scattering should be com- 
parable to those of nuclear scattering. The mosaicity 
which governs secondary extinction should a priori not 
be affected. In order to verify these assumptions the ex- 
tinction correction factor y has been calculated for three 
magnetic reflections using the reflned extinction param- 
eters from the nuclear structure reflnement according to 
the anisotropic FullProf model: 



y 



4sin(2l9)(sin(9/A)2 



(4) 



The calculated values have been compared with the ob- 
served ones, which can easily be deduced from the in- 
tensity ratios at H={) T and i^=0.25 T. The results are 
listed in Tab. [Ill It can be seen that the calculated extinc- 
tion factors are to a greater or lesser extent comparable 
with the observed ones. Nevertheless, the extinction of 
magnetic scattering seems to be underestimated. 



TABLE II: Observed and calculated extinction correction pa- 
rameters for three low-angle magnetic reflections. 



(hkl) 


sin6'/A (A-^) 


FM,ohs (10-^2 cm) 


Vobs 


Veal 


(021) 


0.10595 


4.87(1) 


0.39(7) 


0.46(3) 


(002) 


0.12064 


3.39(1) 


0.47(6) 


0.64(3) 


(023) 


0.20083 


4.72(1) 


0.47(6) 


0.60(3) 



B. Real space magnetization density 

The nuclear structure factors, which have been de- 
duced from the unpolarized neutron experiments, have 
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FIG. 3: (Color online) Observed (squares) and calculated 
(dots) flipping ratios as a function of sinO/X. 



been used to derive the magnetic structure factors from 
the observed flipping ratios by solving Eq.[T]with respect 
to Fm' The individual observed and calculated flipping 
ratios are shown in Fig. [3] as a function of sinO/X. As 
the crystal structure is centrosymmetric the experimen- 
tal magnetization density can directly be reconstructed 
by a Fourier synthesis. Fig. [IJb) shows the projection of 
the observed magnetization density onto the b-c plane, 
while Fig. SJ^a) depicts the unit cell viewed along the a 
axis in order to assign the density peaks to the respective 
atoms. As it has been assumed in the previous section, 
the Coc moments do not get saturated, but significant 
magnetization density is present on V and O sites. While 
the density is quite localized for the V, 01 and 02 sites, 
rather diffuse density can be observed around the 03 site. 
The split density peaks of 01 result from the fact, that 
actually two 01 ions are visible in the projection. Simi- 
larly, the density around the Cog ions seems to be much 
higher compared to the Coc ions, which is due to the fact, 
that two Cos ions are contained in the projection, while 
only one Coc ion is projected. Besides the superexchange 
pathways Cos-02-Coc and Cos-03-Coc an interlayer ex- 
change becomes evident with the non-zero magnetization 
density of V and 01. Further ab initio solid state com- 
putations are currently being performed to simulate the 
spin density map of C03V2O8 and to elucidate the com- 
posite mechanisms of the induced magnetic moments on 
the different O and V sites. This will be the subject of a 
forthcoming publication. — 
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FIG. 4: (Color online) (a) Crystal structure viewed along 
the a axis, (b) Experimental and (c) calculated magnetiza- 
tion density as a projection onto the b-c plane. Contour lines 
defining positive values are drawn as solid lines in 0.05 /is/A^ 
intervals between /isA^ and 0.15 /is/A^ and in 0.4 /i^/A^ 
intervals above. Negative isodensities are represented by bro- 
ken lines in 0.1 /is/A^ steps. 



C. Momentum space density 

The MCPs were extracted by taking the difference 
of the scattered intensities /+ and /~ of the respec- 
tive charge Compton profiles. Due to the fact that the 
more intense charge Compton profiles still exhibit rela- 
tively large values at the outermost measured positions 
of Pz = ±10 a.u. (atomic units), the actual number of 



electrons between -10 a.u. and +10 a.u. is evaluated 
from the profiles interpolated using tabulated data for 
the elements resulting from Hartree-Fock calculations.^^- 
Before summing up the magnetic intensity of each detec- 
tor cell, the data have been corrected for the detector cell 
efficiency, sample absorption and scattering cross-section 
according to Ref. [l^. Furthermore, the energy scale of 
each detector cell has been calibrated by measuring a ra- 
dioactive sample with well known emission energies. 
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The experimental MCPs were folded at Pz=0 to increase 
statistical accuracy by taking the average of each branch. 
The area under each profile has been normalized to the 
number of magnetic electrons per formula unit. With the 
use of iron standards the induced ferromagnetic compo- 
nent can be deduced from the magnetic effect, which is 
the relative contribution of the MCP to the total Comp- 
ton profile 

Mo = • 100%. (5) 

The magnetic effects of the respective directional MCPs 
are listed in Tab. IIIII with their corresponding ferromag- 
netic components induced parallel to the scattering vec- 
tor. In the Py-Pz plane the spin moments monotonically 
decrease with increasing angle between the magnetic field 
direction and the [001] direction, which directly reveals 
the magnetic anisotropy with b being the hard axis of this 
systemj^Si^ Effective beam path length dependent multi- 
ple scattering effects due to the rotation of the sample are 
estimated to be less than a few percent of the spin values 
given in Tab. Hill Fig. [5] shows the normalized observed 

TABLE III: Magnetic effects of the respective directional 
MCPs with a magnetic field H applied along the scattering 
vector. i9 denotes the angle between a MCP and the [001] 
direction. 



MCP 




H{T) 


Mo (%) 


S (/is) 


[001] 





2 


0.541 


0.616 


[023] 


17 


2 


0.488 


0.556 


[012] 


34.7 


2 


0.415 


0.472 


[Oil] 


54.1 


2 


0.274 


0.312 


[032] 


64.3 


2 


0.229 


0.261 


[010] 


90 


2 


0.095 


0.108 


[100] 


90 


0.25 


0.251 


0.287 



MCPs, which reveal similar shapes for the seven inves- 
tigated crystallographic directions. Using all profiles ex- 
cept the one along the [100] direction the two-dimensional 
momentum spin density in the Py-Pz plane has been re- 
constructed by the direct Fourier-transform method j^^i^^ 
The calculation has been performed on a grid with a 
distance of 0.1 a.u. between each point. The result is 
shown as a two-dimensional contour plot in Fig. [6ja). 
Low spin density can be recognized inside the first Bril- 
louin zone (BZ), which extends beyond its border along 
the (010) and (001) directions. In the vicinity of the first 
BZ border the density increases more rapidly with in- 
creasing momentum along (021). Peaks are present at 
(p^,j9^) = (0.35,1.85) and (1.4,0.55). 

D. Correlated refinement in both spaces 

The idea behind correlating the density distributions 
in real and momentum space is that the population of 
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FIG. 5: (Color online) Observed (dots) and calculated (solid 
lines) normalized directional MCPs (shifted vertically in order 
to improve clarity, horizontal lines serve as a guide for the 
eye) . The abscissa pz is taken to be parallel to the respective 
scattering vector. i9 denotes the angle between a respective 
MCP and the [001] direction. 



each spin polarized orbital must be represented likewise 
in the observed densities in both spaces. The fact that 
the MCS method only samples the spin part and the 
PND method both the spin and the orbital part of the 
magnetic moment has been handled in the following way. 
As the observed MCPs have been normalized, the area 
under each profile corresponding to the size of the spin 
moment is not refined. The refined parameters were the 
populations of each spin polarized orbital (in correlation 
with the real space quantities), thus only the shape of 
each profile is refined. The refined magnetic moments 
stem solely from the PND data. In order to analyze the 
observed MCPs theoretical ones have been calculated by 
projecting the square of the ab initio wf onto the respec- 
tive scattering vector. Thereby the symmetry relations 
between the different cluster density distributions in the 
unit cell have to be taken into account, which yields two 
and four symmetrically inequivalent CocOe and CogOe 
clusters, respectively (Tab.HV|). The point symmetries of 
the CocOe and CogOe clusters are 2/m.. and .2.,^'^ which 
correspond to 2/m.. and .2/m. in momentum space. Due 
to the special symmetry of the CogOe density, the pro- 
jections of the different clusters in momentum space are 
invariant for the principal axes and those in the Py-Pz 
plane. In the case of the CocOe clusters, the projections 
onto non-principal axes in the Py-Pz yield different pro- 
files, which need to be averaged. The projected orbitals 
have been convoluted with a Gaussian function having a 
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(au.) 



FIG. 6: (Color online) Reconstructed experimental (a) and 
calculated (b) spin momentum density in the Py-pz plane. 
Contours are drawn in 0.025 /is/(a.u.)^ intervals. White solid 
lines depict the boundary of the first Brillouin zone. 



full width at half maximum of the instrumental resolu- 
tion. As reported previously^ the fact that the projec- 
tion of each MO in momentum space has a characteristic 
shape makes it possible to refine its population (3k which 
is the contribution to the observed MCP: 



JmagiPz) 



(6) 



Here, Xfc(p) denotes a momentum space MO/wf. Simi- 
larly, the populations Pk of the real space MOs can be 
used to deduce the magnetic form factors fx (q) of the re- 
spective elements X by calculating the Fourier transform 
of the atomic spin density: 



OO OO CX) 



/x(q) 



Pki^k,x (i*) exp(27rzqr)dr, (7) 



TABLE IV: Symmetry relations between the two and four 
inequivalent CocOe and CosOe, respectively. 

cluster Co position symmetry relation to ci/si 



ci 

C2 

si 

S2 
S3 

S4 



(0,0,0) 

(o,i^) 

(1,^,1) 



xyz 
xyz 
xyz 
xyz 
xyz 
xyz 



where V^/c,x defines the real space MO k only includ- 
ing the atomic orbitals ^^^x of element X=Co, O. With 
this procedure the observed flipping ratios can be refined 
based on a simple aspheric magnetic form factor model 
deduced from ab initio wf. For the V ions the analytic 
approximation of the form factor^^ has been used. 
Refining the population parameters for each MCP in- 
dividually yields excellent agreement with the observed 
profiles. But since the refinement process exhibits nu- 
merous local minima with significantly varying results, it 
has been considered more reasonable to include all MCPs 
in the refinement despite the magnetic anisotropy. The 
respective population parameters Pk have been refined si- 
multaneously in both spaces together with the magnetic 
moments of Co, V and O by minimizing the function 



X 



4e 



obs 



i,obs 



n j 



[Jn,obs{Pz,j) Jn,cal{Pz,j)] 
2 

^j,obs 



(8) 



-00 —00 —00 



with i and j defining discrete data points of the PND 
and MCS experiment, respectively, and n referring to 
the respective MCPs. The refinement yields fairly good 
agreement expressed by Rmcs=^'^ and Rpnd=^'^ for 
the respective experiments. The refined total magnetic 
moments along the a axis are 

/i(COe) = 1.54(4) 
/i(Co,) = 2.87(3) 
MV)= 0.41(4) fiB 

/i(01) =0.05(5) liiB 

/i(02) = 0.35(5) /IB 

/i(03) =0.36(5) /iB. 

Summing the magnetic moments of all ions in the unit 
cell weighted by their site multiplicity and dividing by 
the number of Co ions yields an averaged magnetiza- 
tion of 3.45 iib/Go^^ . This value shows excellent agree- 
ment with the macroscopic magnetization for i7 = 2 T 
along the a axis reported in Ref. 0. The resulting rel- 
ative orbital populations are listed in Tab. |Vl The re- 
fined parameters were used to calculate the MCPs, which 
are depicted as solid lines in Fig. [5l Fig. [3 shows that 
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FIG. 7: (Color online) Observed (dots) and calculated (solid 
lines) MCPs along two different directions. The dotted and 
dashed lines show the contribution of the CocOe and the 
CosOe cluster, respectively. 



the line shapes of the two respective contributing parts 
(CocOe and CogOe MOs) are different concerning the 
ratio between the value at the peak and at Pz=0 and 
that they vary with the projection angle, which is impor- 
tant for a meaningful fit. From the calculated MCPs 

TABLE V: Refined orbital occupation parameters of the 
CocOe and CosOe clusters. 



orbital Coc Cos 

da:y 0.27(2) 0.12(2) 

d:,^ 0.27(2) 0.12(2) 

dy^ 0.16(2) 0.26(2) 

d^2_y2 0.17(2) 0.30(2) 

ds,2_,2 0.13(2) 0.20(2) 



the momentum space spin density has been reconstructed 
[Fig.EJb)]. The calculated real space magnetization den- 
sity map [Fig. [11(c)] has been obtained by a Fourier syn- 
thesis of the calculated magnetic structure factors. The 
main features of the respective density maps coincide 
well, although some differences are evident: the dip in the 
momentum space density around Pz=0 is not pronounced 
well in the calculated map and has a shape, which is ro- 
tated by 90° with respect to the observed map. This 



possibly results from strong hybridization effects between 
the Co3d and 02p orbitals. Furthermore, the disagree- 
ments between the experimental and calculated MCPs 
are most probably due to the limitation of the cluster 
calculations. Considering a larger sized cluster or the 
solid state, the gap between the experiment and fitting 
could be reduced. The real space spin density of the 01 
and 03 sites is slightly underestimated. Furthermore, 
density peaks exist next to the spine site density along 
the z axis, which do not coincide with atomic positions. 
However, this fact can be attributed to truncation effects 
in the Fourier series. 



E. Discussion 

The study presented here combining several meth- 
ods reveals very interesting magnetic properties of the 
kagome staircase system C03V2O8. The previous as- 
sumption, that the ferromagnetic structure at zero mag- 
netic field is not fully ordered because of the Coc only 
exhibiting 1.54 /i^, is disproved. Previous macroscopic 
magnetization measurements'^ indeed showed a saturated 
moment of approximately 3.4 /j^b per Co site at H=2 T 
along a, but the results of the polarized neutron single- 
crystal diffraction experiment with adequate extinction 
correction presented here reveal that the field dependent 
increase of magnetization stems from the V, 02 and 03 
sites. Hereby, the V and 02 site show quite localized 
magnetization density, while the 03 density seems to be 
smeared out due to truncation effects in the Fourier se- 
ries. A periodic ab initio calculation confirms the exis- 
tence of magnetization density on the V and O sites and 
will be presented elsewhere. The spin polarized density 
on 02 and 03, which are those oxygen ions in the CocOe 
clusters, may be a strong indication for a partially co- 
valent character of the Coc ions and the reason for their 
relatively low magnetic moment compared to C05. The 
magnetization density distribution clearly exhibits the 
superexchange pathways between the two different Co 
sites, but it indicates also the interlayer coupling, which 
is mediated by the V-01 bridge. Combining the meth- 
ods of polarized neutron diffraction and magnetic Comp- 
ton scattering allowed us to refine the occupations of the 
Co3d orbitals in a stable way. Like it has been previ- 
ously reported but with inverted values^^ the two crys- 
tallographically different Co ions exhibit different spin 
polarized orbital occupations. While the unpaired elec- 
trons are equally distributed between the t2g and Cg levels 
for the Cos ion, the magnetic signal stems by only 30% 
from the Cg orbitals for the Coc ion as a consequence of 
the spin transfer from the surrounding O ions. Concern- 
ing the Cg orbitals of both ions the basal plane orbital 
dx2_y2 is more populated than the apical d^z'^-r'^ orbital. 
This possibly indicates a higher exchange interaction be- 
tween the Cog ions via an intermediate 02 ion. In the 
case of the Coc ions it could be a hint that the magnetic 
exchange with the spine Co ions takes place preferentially 
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via an 03 ion. 
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